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1 Introduction 

In this report we estimate the sensitivity of experiments optimized for measuring 
|/7e3p and the sign of Am^g in various Fermilab main injector beam configurations. 
We calculate the sensitivity to —^ appearance making the assumption of 3 neu- 
trino generations, and that only solar and atmospheric neutrino oscillations signals are 
"real." We have implemented an exact calculation of oscillation probabilities in matter 
from Barger et.al. for the FNAL-Soudan (732km), and hypothetical FNAL-BNL 
(1500km) and FNAL-SLAC (2900km) beamhnes, as shown in figures |l[ |, and |. Of 
course, these appearance probabilities must be folded into the neutrino flux spectrum 
produced by a particular Fermilab main injector neutrino beam configuration, as well 
as the cross-section in the far detector. For the neutrino flux spectrum we have used 
the GEANT simulation of the NuMI beamline |^ , and for the neutrino cross-section 
we have used the Soudan 2 |^ event generator. For the purposes of illustration, 
we shall use the test point Amjg — 0.003ey^ (regular hierarchy), |/7e3p = 0.003 
(ie. = 3°), [7^3 = C/23, Am2^,„, = 0.00003ey2, f/,^ = C/f^, and phase = 0. 
Using these parameters, the numbers of i'^ — > v^, charged current interactions per 
kiloton-year (in steel) are as follows: 



Event rates for steel exposure 

regular hierarchy, Am^g = 0.0036^^, |C/e3p = 0.003 




FNAL-Soudan 


FNAL-BNL 


FNAL-SLAC 


PH2 low v^, -> Ve CC 


0.93 / kt-yr 


0.55 / kt-yr 


0.23 / kt-yr 


PH2 medium v^J_ v,, CC 


1.41 / kt-yr 


1.15 / kt-yr 


0.79 / kt-yr 


PH2 high Z/^ ^ i^e CC 


0.98 / kt-yr 


0.95 / kt-yr 


0.87 / kt-yr 



Note that the peak energy values for the PH2 low, medium, and high energy 
beams are roughly 3, 6, and 12 GeV respectively. In this report we will consider two 
examples of detectors: a 10-kt "MINOS-like" detector, and a 10-kt "OPERA-hke" 
detector optimized for z/g appearance. We define "MINOS-like detector" in this report 
to be a 10-kt detector composed of alternating 2.54cm Fe plates and scintillator strips. 
We define the "OPERA-like" detector as follows. 

In an earlier note, we have described how a hybrid emulsion detector (HED) can 



obtain low background measurements of 



and Vu 



oscillation events 



in the NuMI beam g. We now consider a HED designed exclusively to measure 
Vf,. The tracking requirements for v-r appearance are more stringent than those for 
Ve appearance, and thus we can significantly relax certain parameters of a HED 
optimized for Vt § if we only care about identifying Vf,. The main difference is that 
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we propose to use 5 mm Fe instead of 1 mm Pb as in OPERA Thus, for roughly 
the same emulsion cost as in OPERA (a 2-kt detector) we obtain a detector with 3.5 
times greater mass. 

The detector would be composed of 10 cm thick steel-emulsion stacks, separated 
by RPC planes. The RPC planes would be used for triggering and event location. 
Emulsion sheets would be sandwiched between 5 mm thick steel plates. If the planes 
have an area of 8x8 square meters, then the mass of one steel-emulsion target plane 
would be 50 tons. The entire detector would consist of 200 planes. Note that there is 
enough floor space behind the first two MINOS supermodules in the MINOS Soudan 
cavern for such a detector. 

We now attempt to estimate how much a HED would cost. We use the cost of 
emulsion film from the OPERA proposal (converted to US dollars). We estimate 
the cost of salary, wages, and institutional overhead to be the same fraction as in 
the MINOS TDR [||, ie 50% of total cost before contingency. We add a 50% contin- 
gency. This results in a total detector cost of roughly $18. 4M per kt. The rough cost 
breakdown is shown in the following table: 



Cost breakdown for 10-kt HED 


component 


cost per unit 


^ units 


total cost 


emulsion film 


$150 / 


256,000 


$38.4M 


steel 


$1000 / ton 


10,000 


$10.0M 


RPC-|-electronics 


$500 / 


25,600 


$12. 8M 


wages, salary, overhead 






$61. 2M 


contingency 






$61.2M 


total 






$183. 6M 



We now turn to a comparison of MINOS-like and OPERA-like detector analyses. 

2 MINOS-like data analysis 

The MINOS-like detector analysis (taken from 0) is based upon the full "official" 
GEANT simulation of the MINOS detector g. We use the case of PH2 medium 
to SLAG as an example. The basic strategy can be summarized as follows: 1) re- 
ject events with > IGeV, 2) fraction of energy in the highest energy cluster 
EcLUST / Etot > 0.7 (see figure ||), 3) number of strips in the highest energy cluster 
N ST RIPS >— d (see figure ^), 4) neural net estimator consistent with i/g (see fig- 
ure H), and 5) reject events in which the total energy {Etot) does not fall within 
some energy range optimized to Am^ (see figures ^ and |^) . The numbers of events 
after each sequential cut are tabulated in table |^. 

3 HED data analysis 

The HED data analysis Q is based upon simple gaussian smearing of individual 
Monte Garlo particle truth information. Of course, the smearing is appropriate to 
the HED detector under consideration. The basis for the inputs to the smearing 
come from the OPERA proposal (which has an ^0=0. 18 for each target plate), 
with the appropriate extrapolation to the thicker plates in the HED described in this 
report (Aro=0.28). The resolution on electromagnetic showers from track counting in 
the OPERA proposal is 20%/sqrt(E), so we use a gaussian smearing of 25%/sqrt(E) 
for electrons and gammas in our calculation. The resolution on charged particle 
momentum using the multiple Goulomb scattering method in the OPERA proposal is 
16%, so we use a gaussian smearing of 20% to (non-electron) charged particles in our 
calculation. We use a 1 milliradian angular resolution for charged particle direction 
(assuming that we use the emulsion film described in the OPERA proposal) and we 
add in quadrature the angular smearing due to multiple coulomb scattering in the 
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target steel plate (using a randomly chosen scattering vertex position in the target 
plate). 

We use the case of PH2 medium to SLAC as an example. The data reduction 
can be summarized as follows: 1) reject events with > IGeV, 2) require the 
highest energy electromagnetic shower (either from an electron or prompt gamma) 
to be greater than some optimized threshold (sec figure^), 3) reject prompt gammas 
which do not convert before passing through the first emulsion sheet after the primary 
interaction vertex (this rejects 88% of prompt gammas), 4) require missing transverse 
momentum^ to be less than some optimized threshold (see figure |l^), and 5) require 
the total visible energy to fall within an energy range optimized to Am^ (see figures |ll| 
and |l^) . The numbers of events after each sequential cut are tabulated in table 

4 Results 

We assume 40 kt-yr exposures of the detectors in Fermilab main injector beams whose 
fluxes have been upgraded by a factor of 4. We also assume that the systematic error 
is dominated by uncertainty in the number of background events, and that this error 
is 10%. With these assumptions, we may summarize the significance of the signal in 
units of astat and agyst after all data reduction cuts in the following table: 



Significance of the signal after all cuts in units of astat {<^syst) 
regular hierarchy, Am^g = 0.0036^^^ |C/e3p = 0.003 




MINC 
Soudan 


)S-like de 
BNL 


;ector 
SLAC 


OPE] 

Soudan 


lA-like dc 
BNL 


itector 
SLAC 


PH2 low 


1.5 (0.8) 


2.3 (2.1) 


1.4 (2.2) 


2.6 (2.2) 


3.2 (4.3) 


2.0 (5.2) 


PH2 medium 


1.1 (0.4) 


1.9 (1.3) 


3.1 (2.6) 


1.4 (1.3) 


2.5 (4.2) 


4.5 (10.0) 


PH2 high 


0.4 (0.1) 


0.7 (0.4) 


1.7 (0.8) 


1.1 (0.8) 


1.9 (2.6) 


3.3 (7.5) 



The predicted sensitivities to iC/esI^ as a function of Atojs are shown in figures 
through |l^. Two general conclusions may be made: 1) OPERA-like is better than 
MINOS-like, and 2) running "on the highest energy oscillation peak" is best. If 
we take these calculations at face value, then the best beam/site/detector setup for 
Ato23 = 0.003eT^^ appears to be the PH2 medium beam pointed at SLAC/HED. 
(If Am^g is slightly higher or lower, then a slightly higher or lower energy beam 
is optimal.) The predicted sensitivity to |C/e3p including statistical and systematic 
errors for a 40-kt exposure of HED at SLAC in PII2 medium (4x regular flux) is 
shown in figure Similarly, the discovery potential for |C/e3p and sign(Am23) is 
shown in figure |2C| . Note that appearance signals and the sign of Am^g can 

be obtained at the 3-sigma level for 6*13 >3 degrees. 

We would like to thank Stanley Wojcicki, Vittorio Paolone, and Deborah Harris for 
several useful comments and discussions. Also, we would like to thank Mark Messier 
for providing a cross-check of our implementation of the Barger et.al. method for 
calculating neutrino oscillations in matter. 

This work was sponsored by the National Science Foundation. 

References 

[1] Matter effects on three-neutrino oscillations. V. Barger, et.al. 
Phys.Rev.D22, pp.2718-2726 (1980). 

[2] Conceptual design for the techincal components of the neutrino beam for 
the main injector (NuMI). J.Hylen et.al. Fermilab-TM-2018. (Sept. 1997) 

^Note that in order to minimize the missing transverse momentum smearing due to multiple 
Coulomb scattering we only use charged particles with momentum greater than 400 MeV in the two 
independent transverse momentum sums. 



3 



[3] New results on atmospheric neutrinos from Soudan 2. W.A.Mann for 
Soudan-2 Collaboration. PDK-752, TUHEP-00-05, (Jul 2000). Presented 
at 19th International Conference on Neutrino Physics and Astrophysics 



Neutrino 2000, Sudbury, Ontario, Canada, 16-21 Jun 2000. hep- 



ex/0007031 



The MINOS Detectors Technical Design Report. MINOS CoUaboration. 
NuMI-L-337. 

Optimization of a hybrid emulsion detector for MINOS. B.Bcrnevig et.al. 
NuMI-E-525 (Aug. 1999) 

Calculation of v^. neutrino oscillation signals in a MINOS hybrid 

emulsion detector. L.Wai and S.Wojcicki. NuMI-E-547 (Oct. 1999) 

Study of Vf, -> sensitivity in MINOS. L.Wai et.al. NuMI-L-714. 
(Dec.2000) 

MINOS detector GEANT simulation. R.Hatcher. See rhatcher@fnal.gov. 

An appearance experiment to search for — > Vt oscillations in the CNGS 
beam. OPERA collaboration. CERN/SPSC 2000-028 (July 2000) 

Limits on Neutrino Oscillations from the CHOOZ Experiment. CHOOZ 
Collaboration. (19 July, 1999) |hep-ex/990703^ 

Physics at a Neutrino Factory. FERMILAB-FN-692 (April 2000) 



4 





Figure 1: Probability for v^, -> Ve for FNAL-Soudan (L=732km), given |(7e3p = 0.003 
(ie. 6»i3 = 3°), Ato^3 = 0.003ey^ Ato?2 = 0.00003ey^ U^^^ = U^^, and U^^ = U^^- 
The black lines are for neutrinos, and the red lines are for anti-neutrinos. The panels 
on the left side have regular mass hierarchy, and the panels on the right side have 
inverted mass hierarchy. The upper panels have zero phase, the middle panels have 
+90 degrees phase, and the lower panels have -90 degrees phase. 
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Figure 2: Probability for i/^ for FNAL-BNL (L=1500km), given \Ues.\^ = 0.003 

(ie. 6»i3 = 3°), Ato^3 = 0.003ey^ Amfa = 0.00003ey2^ L^^^ = U^^, and C/^^ = U^^- 
The black lines are for neutrinos, and the red lines are for anti-neutrinos. The panels 
on the left side have regular mass hierarchy, and the panels on the right side have 
inverted mass hierarchy. The upper panels have zero phase, the middle panels have 
+90 degrees phase, and the lower panels have -90 degrees phase. 
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Figure 3: Probability for ly^ for FNAL-SLAC (L=2900km), given |C/e3p = 0.003 
(ie. 6li3 = 3°), Ato^3 = 0.003ey^ Ato?2 = 0.000036^^^ = U^^, and C/^^ = U^^- 
The black lines are for neutrinos, and the red lines are for anti-neutrinos. The panels 
on the left side have regular mass hierarchy, and the panels on the right side have 
inverted mass hierarchy. The upper panels have zero phase, the middle panels have 
+90 degrees phase, and the lower panels have -90 degrees phase. 
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FNAL-SLAC (2900 km) 40 kt-yr exposure of MINOS-like detector in PH2 medium (4x) 
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CC 
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NC 


sequential cuts: 




beam Ve 


Vfi Vpi 


Vfi Vt 


< lOGeV 


E^ > lOGeV 


all events 
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4589.8 


1325.7 


3312.7 
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FNAL-SLAC (2900 km) 40 kt-yr exposure of HED detector 


in PH2 medium (4x) 
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Figure 4: PH2 mcdium/MINOS-likc/SLAC event distributions of the fraction of total 
energy in the highest energy cluster {Eclust / Etot)- We require Eclust/Etot > 
0.7. The histograms are normalized to 4x 40 kt-yr data samples. An exact calculation 
of the oscillation probability in matter was used for Am^s = 0.0036^^ (regular mass 
hierarchy), |[/e3p = 0.003 (ie. 6^13 = 3°), U^^ = U^^, U^^ = U^^, Amf^ = 0.00003ey2, 
and phase ^ = 0. 
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Figure 5: PH2 mcdium/MINOS-likc/SLAC event distributions for number of strips 
(NsTRiPs) in the highest energy cluster. We require Nstrips >= 9. The his- 
tograms are normaUzed to 4x 40 kt-yr data samples. An exact calculation of the 
oscillation probability in matter was used for Am^g = 0.003eV^ (regular mass hierar- 
chy), |[/e3p = 0.003 (ie. 013 = 3°), Uf^^ = U^^, U^^ = U^^, Amj^ = 0.00003ey2, and 
phase (f) = 0. 
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Figure 6: PH2 medium/MINOS-like/SLAC event distributions for the neural net 
estimator (Y). We require y > 0. The histograms arc normahzcd to 4x 40 kt-yr data 
samples. An exact calculation of the oscillation probability in matter was used for 
Am^3 = O.OOSey^ (regular mass hierarchy), |C/e3|^ = 0.003 (ie. 6»i3 = 3°), U^^ = U%, 
Uei = Ue2, Amfa = 0.00003ey2, and phase <?!> = 0. 
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Figure 7: PH2 mcdiuni/MINOS-like/SLAC event distributions for the total energy 
{Etot)- We require 400 photoelectrons < Etot < 800 photoelectons. The his- 
tograms are normalized to 4x 40 kt-yr data samples. An exact calculation of the 
oscillation probability in matter was used for Am^g = O.OOSey^ (regular mass hierar- 
chy), |[/e3p = 0.003 (ie. 013 = 3°), Uf^^ = U^^, U^^ = U% Amj^ = 0.00003ey2, and 
phase (f) = 0. 
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Figure 8: PH2 medium/MINOS-like/SLAC total energy event distributions for back- 
ground (black line) and background plus signal (red line). We require the total energy 
to be in the range between 400 and 800 photoelectons (green region) . The histograms 
are normalized to 4x 40 kt-yr data samples. An exact calculation of the oscilla- 
tion probability in matter was used for Am^g = 0.003eV^ (regular mass hierarchy), 
|C/e3p = 0.003 (ie. 9i3 = 3°), U^^ = U^^, U^^ = U^^, Amj^ = O.OOOOSeT^^^ and phase 
(/) = 0. 
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Figure 9: PH2 medium/HED/SLAC event distributions for the maximum energy 
electromagnetic shower (Ee^-y). We have used a 25%/sqrt(E) electromagnetic shower 
energy resolution (from track coimting). Wc require -Ee,7 > 5 GcV. The histograms 
are normalized to 4x 40 kt-yr data samples. An exact calculation of the oscilla- 
tion probability in matter was used for Am^^ = O.OOSeV'^ (regular mass hierarchy), 
|C/e3p = 0.003 (ie. 6»i3 = 3°), U^^ = U^^, U^^ = U%, Am?2 = O.OOOOSeF^, and phase 
^ = 0. 
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Figure 10: PH2 medium/HED/SLAC event distributions for the missing transverse 
momentum (PT) . We have used a 20% track momentum resolution (from the muhiple 
scattering method) and 1 mrad track angular resolution. (Obviously, we have not used 
neutrinos, neutrons, or K-longs in the Px,Py sums.) We have incorporated multiple 
coulomb scattering; in order to minimizing the smearing due to multiple scattering 
we only use charged particles with momentum greater than 400 MeV in the sums. We 
require missing transverse momentum PT < 400 MeV. The histograms are normalized 
to 4x 40 kt-yr data samples. An exact calculation of the oscillation probability in 
matter was used for Am^g = O.OOSeT^-^ (regular mass hierarchy), |J7e3p = 0.003 (ie. 
6»i3 = 3°), = U%, = U% ^ml2 = 0.00003ey2, and phase = 0. 
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Figure 11: PH2 medium/HED/SLAC event distributions for the visible energy {Eyis). 
We have used a 20% track momentum resolution (from the multiple scattering 
method) and a 25%/sqrt(E) electromagnetic shower energy resolution (from track 
counting). We require Eyis < 10 GeV. The histograms are normalized to 4x 40 kt- 
yr data samples. An exact calculation of the oscillation probability in matter was 
used for Am^g = O.OOSeV^ (regular mass hierarchy), jfJeal^ = 0.003 (ie. 613 = 3°), 
Uu3 = U^3, U^i = U^2, Amfa = O.OOOOSeF^, and phase = 0. 
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Figure 12: PH2 medium/HED/SLAC total energy event distributions for background 
(black line) and background plus signal (red line). We require the total energy to be 
less than 10 GeV (green region). The histograms are normalized to 4x 40 kt-yr data 
samples. An exact calculation of the oscillation probability in matter was used for 
Amis = 0.003eT/2 (regular mass hierarchy), |f7e3p = 0.003 (ie. Oi^ = 3°), C/^^ = u^^, 
Uei = U^2, Am?2 = 0.00003ey2^ and phase (j) = 0. 
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Figure 13: (4x) 40 kt-yr 90% confidence level limits on |J7e3p with Fermilab main 
injector beam pointed at Soudan/MINOS-like. We have assumed a systematic error 
of 10% on the number of background events. An exact calculation of the oscillation 
probability in matter was used for regular mass hierarchy, U^^ — U^^, U^i 
Ato?2 = 0.00003ey2^ and phase = 0. 
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90% C.L. limits, regular hierarchy, U ^3=U^3 
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Figure 14: (4x) 40 kt-yr 90% confidence level limits on |?7e3p with Fermilab main 
injector beam pointed at BNL/MINOS-like. We have assumed a systematic error of 
10% on the number of background events. An exact calculation of the oscillation 
probability in matter was used for regular mass hierarchy, U^^ — C/^g, U^-^ 
Ato?2 = 0.00003ey2^ and phase = 0. 
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Figure 15: (4x) 40 kt-yr 90% confidence level limits on |J7e3p with Fermilab main 
injector beam pointed at SLAC/MINOS-like. We have assumed a systematic error 
of 10% on the number of background events. An exact calculation of the oscillation 
probability in matter was used for regular mass hierarchy, t/^g — U^^, U^^ 
and phase (j) — 0. 
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Figure 16: (4x) 40 kt-yr 90% confidence level limits on |?7e3p with Fermilab main 
injector beam pointed at Soudan/HED. We have assumed a systematic error of 10% on 
the number of background events. An exact calculation of the oscillation probability 
in matter was used for regular mass hierarchy, C/^g = U^^, U^-^ = U^2^ ^"^12 = 
O.OOOOSey^ and phase = 0. 
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Figure 17: (4x) 40 kt-yr 90% confidence level limits on |?7e3p with Fcrmilab main 
injector beam pointed at BNL/HED. We have assumed a systematic error of 10% on 
the number of background events. An exact calculation of the oscillation probability 
in matter was used for regular mass hierarchy, C/^g = U^^, U^-^ = U^2^ ^"^12 = 
O.OOOOSey^ and phase = 0. 



23 



SLAC/HED sensitivities! 



> 



CM CM 



E 0.005 

< 



0.004 



0.003 



0.002 — 



0.001 



90% C.L. limits, regular hierarchy, U i^3=U^ 



(i3 t3 



PH2 low SLAC/HED 
PH2 med SLAC/HED 
PH2 high SLAC/HED 




logiodUesI) 



Figure 18: (4x) 40 kt-yr 90% confidence level limits on |?7e3p with Fermilab main 
injector beam pointed at SLAC/HED. We have assumed a systematic error of 10% on 
the number of background events. An exact calculation of the oscillation probability 
in matter was used for regular mass hierarchy, C/^g 
O.OOOOSey^ and phase = 0. 
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Figure 19: Predicted 90% confidence level limits on |?7e3p with Fermilab main injector 
beam. We have assumed a systematic error of 10% on the number of background 
events. An exact calculation of the oscillation probability in matter was used for 
regular mass hierarchy, f/^g = U^^, U^^ = U^2^ ^'^12 — O.OOOQSeT^^, and phase 
(h = 0. 
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Figure 20: Discovery potential for [C/esl^ ^^id sign(Am23) with Fermilab main injector 
beam. We have assumed a systematic error of 10% on the number of background 
events. For sign(Am23) we have assumed an additional sample of data with the same 
exposure except with an anti-neutrino beam. An exact calculation of the oscillation 



probability in matter was used for regular mass hierarchy, C/^g 



Am?2 = 0.00003ey^ and phase 
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